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^ ! Abstract 

We study acoustiophonon-induced relaxation of charge excitations in single and tunnel-coupled 
quantum dots containing few confined interacting electrons. The Full Configuration Interaction 
approach is used to account for the electron-electron repulsion. Electron-phonon interaction is 
accounted for through both deformation potential and piezoelectric field mechanisms. We show 
that electronic correlations generally reduce intradot and interdot transition rates with respect to 
corresponding single-electron transitions, but this effect is lessened by external magnetic fields. On 
the other hand, piezoelectric field scattering is found to become the dominant relaxation mechanism 
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phonon coupling in properly designed single-electron quantum dots are shown to hold also in 
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as the number of confined electrons increases. Previous proposals to strongly suppress electron- 



multi-electron devices. Our results indicate that few-electron orbital degrees of freedom are more 
stable than single-electron ones. 
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I. INTRODUCTION 



Recent experiments have appointed coupling to acoustic phonons as the main source of 
electron relaxation in the excited states of semiconductor quantum dots (QDs) with weak 
lateral confinement^ and coupled quantum dots (CQDs) with small tunneling energy,—^ 
i.e., with the electronic relevant energy scale in the few-meV range. This usually restricts 
lifetimes to the order of nanoseconds,— imposing severe limitations on the performance of 
QD devices which rely on the coherent dynamics of electron charge states. A prominent 
example are charge qubits, whose main disadvantage, as compared to spin qubits, is the 
faster decoherence rate. 1 

As a first step to eventually gain control over charge relaxation rates, several theoret- 
ical works have investigated the fundamental physics of electron-phonon coupling in QD 
structures. Bockelmann described the qualitative effect of lateral (spatial and magnetic) 
confinement on the electron transition rates in QDs.— The domains of deformation potential 
(DP) and piezoelectric field (PZ) interactions, the two main scattering mechanisms leading to 
interaction with acoustic phonons, were established for both QD and CQD structures . 7 ! 8 ' 9 ' 10 
Building on these works, methods to minimize electron-phonon coupling in QDs and ver- 
tically CQDs were recently proposed which may bring electron lifetimes in the range of 
microseconds, thus making them comparable to usual spin relaxation rates . 7 ' 11 ' 12 ' 13 

Noteworthy, all of the aforementioned experimental and theoretical works studied charge 
relaxation at a single-electron (SE) level. So far, multi-electron (ME) systems have been 
essentially restricted to the context of spin relaxation , 1 ' 2 ' 14 ' 15 and only lately charge relax- 
ation processes in ME systems started to be considered . 16 ' 17 In particular, in Ref. ^ we 
have investigated Coulomb correlated systems, and we have reported numerical evidence 
that electronic correlations generally lead to reduced decay rates of the excited electronic 
states in weakly confined dots, as compared to independent-particle estimates. This finding 
suggests that ME devices might better exploit the discrete energy spectrum of QDs. 

In this paper, we extend our previous investigation (Ref. Il 71 ) analyzing in detail the 
mechanisms by which Coulomb interaction affects charge relaxation rates in ME QD systems. 
We also investigate the effect of external magnetic fields on intradot transition rates, and 
isospin transition rates in ME vertically CQDs. Electron-phonon interaction is included 
through both DP and PZ scattering channels (the latter was neglected in Ref. |l7| as well as 
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in many theoretical investigations). 

The paper is organized as follows. In Section |Il] we describe our theoretical model. In Sec- 
tion lllll we study ME charge relaxation in single QDs as a function of the lateral (spatial and 
magnetic) confinement, and interpret the correlation- induced scattering reduction in terms 
of the SE configuration mixing. The applicability in ME QDs of magnetic-field-based sug- 
gestions for controlling charge relaxation rates, previously reported for SE structures ' 12 ^ 3 , 
is addressed in this section. In Section |IV] we investigate the effect of interdot distance on 
the isospin transition rate of ME vertically CQDs, and finally in Section IS we summarize 
our results. 

II. THEORETICAL FORMULATION 

The theoretical model we use is similar to that described in Ref. Q for single-electron 
excitations, but now considering iV-electron states (N = 1 — 5). We study disk-shaped QDs 
where the lateral confinement is much weaker than the vertical one, and the dot and sur- 
rounding barrier are made of materials with similar elastic properties. Several QD structures 
reported in the literature fit this description . 1 i 2 i 18 

A number of useful approximations can be made for such QDs. First, since the weak 
lateral confinement gives inter-level spacings within the range of few meV, only acoustic 
phonons have significant interaction with free carriers, while optical phonons can be safely 
neglected. Second, the elastically homogeneous materials are not expected to induce phonon 
confinement, which allow us to consider only bulk phonons. Finally, the different energy 
scales of vertical and lateral electronic confinement allow us to decouple vertical and lateral 
motion in the building of SE spin-orbitals. We then take a parabolic confinement profile 
in the in-plane (x,y) direction, with single-particle energy gaps Hlvq, which yields the Fock- 
Darwin states.-^ In the vertical direction (z) the confinement is provided by a rectangular 
quantum well of length L z , with the barrier height determined by the band-offset between 
the QD and barrier materials. The quantum well solutions are derived numerically. Spin- 
orbit coupling is neglected in this work, since the long spin relaxation times measured in 
QD structures similar to those we study here indicate that the spin degrees of freedom are 
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well separated from the orbital ones. 1 ' Therefore, the SE spin-orbitals can be written as: 



where 4> n ,m is the n— th Fock-Darwin orbital with azimuthal angular momentum m, £ g is 
the symmetric quantum well solution with parity g with respect to the reflection about the 
z = plane, with g = (1) denoting even (odd) parity, and \ a is the spinor eigenvector of 
the spin z— component with eigenvalue o [a = ±1/2). We will also label Fock-Darwin states 
with the standard notation nl, where I = s, p±, d± . . . correspond to m — 0, ±1, ±2 . . ., 
respectively. 

As for the inclusion of Coulomb interactions, we need to go beyond the Hartree-Fock ap- 
proximation in order to include electronic correlations, which are critical in phonon-induced 
electron scattering processes.— Moreover, since we are interested in the relaxation time of 
excited states, we need to know both ground and excited states with comparable accuracy. 
Our method of choice is the Full Configuration Interaction approach: the ME wave func- 
tions are written as linear combinations |\I/ a ) = Yli c <n\&i) > where the Slater determinants 
= are obtained by filling in the SE spin-orbitals a with the N electrons in 

all possible ways consistent with symmetry requirements; here creates an electron in the 
level a. The ME ground and excited states can then be labeled by the total angular mo- 
mentum z component M = 0, ±1, ±2 . . ., total parity G = 0,1, total spin S, and total spin 
z-projection S z . The fully interacting Hamiltonian is numerically diagonalized, exploiting 
orbital and spin symmetries.— 

We assume zero temperature, which suffices to capture the main features of one-phonon 
processes.— At this temperature, only phonon emission processes are present. We evaluate 
the relaxation rate between selected initial (occupied) and final (unoccupied) ME states, b 
and a, using the Fermi golden rule: 



where the electron states \^f a ) i a — a ,b) have been written explicitly as linear combina- 
tions of Slater determinants, V vai is the interaction operator of an electron with an acoustic 
phonon of momentum q via deformation potential [v = DP) or piezoelectric field [y = PZ) 
interaction, E a stands for the a electron state energy and E q represents the phonon energy. 



if) a (x, y, z; a) = <f) n>m (x, y) i 9 {z) X, 



(1) 
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The electron-phonon interaction matrix element can be written more explicitely as 



($,|K q |$,> = MM) {H e^ qr |^>, (3) 

where the right-most term is the electron form factor and M w (q) is a prefactor which depends 
on the scattering mechanism v. 7 It is worth noting that, whereas for DP scattering Mp,p oc 
\/iqJ, for PZ scattering Mpz oc 1/yfqj. As a result, DP scattering is dominant when 
the emitted phonon energy is sufficiently large, while PZ scattering dominates at small 
phonon energy. For SE transitions in GaAs QDs, even in the weakly-confined regime, the 
DP mechanism usually prevails in the absence of external fields. However, we have recently 
shown 7 that the PZ mechanism may rapidly become dominant in the presence of a vertical 
magnetic field which tends to suppress the single-particle gaps. 

One can see from the above expressions that Coulomb interaction influences electron 
scattering with phonons in two ways. First, it introduces changes in the electron energy 
gaps Eb — E a , and hence in the energy and momentum of the emitted phonon.— Second, 
it introduces changes in the orbital part of the electron state, hence changing the electron- 
phonon wave functions coupling. The latter effect is reflected in Equation (J2J) through the 
Slater determinant coefficients. Indeed, since the total scattering rate is but a weighted 
sum of SE contributions, the general behavior of SE scattering events will also apply to the 
ME case. However, the weight of each SE contribution strongly depends on the number of 
particles, the regime of correlations, and the presence of external fields, so that important 
changes in the ME relaxation rates should be expected when varying these parameters. 

In this work we consider mostly relaxation rates corresponding to the fundamental spin- 
conserving transition in single and coupled QDs, i.e. transitions involving the ground state 
and the first excited state with the same (S, S z ) quantum numbers. This transition could 
be monitored, e.g., by means of pump-and-probe techniques,— "2»i4 since relaxation to or from 
any intermediate state with different spin should be much slower and therefore it will barely 
interfere. 

Below we shall investigate GaAs/Alo.3Gao.7As QDs, using the following material 
parameters?^ electron effective mass m* = 0.067, band-offset V c = 243 meV, crystal density 
d = 5310 kg/m 3 , acoustic deformation potential constant D = 8.6 eV, effective dielectric 
constant e = 12.9, and piezoelectric constant hu = 1.41 • 10 9 V/m. For the sound speed 
C/j , we take into account that in cylindrical QDs most of the scattering arises from phonon 
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propagation close to the growth direction.— We then assume that the QDs are grown along 
the [1 0] direction and use the corresponding values c LA = 4.72 • 10 3 m/ s and c TA = 3.34- 10 s 
m/s.— In our calculations we deal with QDs with lateral confinement energies which in some 
cases are rather weak (Huq < 1 meV). Well-converged few-body states are obtained for such 
structures using a basis set composed by the Slater determinants which result from all pos- 
sible combinations of 62 SE spin-orbitals with N electrons. Due to the strong confinement 
in the vertical direction, only the lowest g = (for single QDs) or the lowest g — 0, 1 (for 
CQDs) eigenstates are included in the single-particle basis. 

III. SINGLE QUANTUM DOTS AND MAGNETIC FIELD 

In this section we first study charge relaxation rates in QDs with iV interacting electrons 
as a function of the harmonic lateral confinement originated by electrostatic fields, and next 
consider the effect of adding a magnetic field. To study the lateral confinement, we vary the 
characteristic frequency of the confining parabola in the huo ~ — 6 meV range, thus moving 
from a strongly- to a rather weakly-correlated regime. Figure ^ depicts the corresponding 
results for N = 1 — 5 electrons in QDs with height L z = 10 nm. It can be observed that for 
all N the qualitative shape of the relaxation rate curve is similar to that of the SE case (top 
panel): it shows two maxima, connected with the PZ and DP scattering mechanisms, and it 
vanishes at small and large confinement energies due to the small phonon density and small 
electron-phonon coupling, respectively— Another trend observed in Fig. ^ is the shift 
of the scattering rate maxima towards larger confinement energies with increasing number 
of electrons, as well as the increasing relative height of the PZ maximum. Both features 
follow from the increasing density of states with larger N, which leads to smaller inter-level 
spacing and, therefore, to larger hujQ values associated with a given phonon energy. We 
illustrate this in Fig. Q with downward arrows to point at the confinement energies which 
give two selected phonon energies, E q = 1.3 (solid arrowhead) and E q = 2.0 meV (empty 
arrowhead). It can be seen that there is a shift towards larger Tlujq values with increasing 
iV.— An important implication of this result is that PZ scattering, which is negligible except 
at very weak confinement energies in the SE picture, may actually become the dominant 
scattering mechanism for usual QD confinement energies in the ME case (see, e.g., hu>o = 2 
meV in the N = 5 picture). In light of this, some of the estimates in previous investigations, 
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which studied electron-phonon coupling in weakly- confined QDs considering DP interaction 
only, may need a revision . 16 1 17 

One also observes in Fig. Q that the phonon energy which gives maximum scattering 
is approximately constant, regardless of the number of electrons and lateral confinement. 
For example, E q pa 1.3 meV for the DP maximum. This indicates that the scattering is 
mainly determined by the electron-phonon coupling along the vertical direction: indeed, at 
E q pa 1.3 meV the longitudinal acoustic phonon wavelength gives maximum coupling with 
the electron wave function in the quantum well, which does not depend on Huq and is weakly 
affected by Coulomb interactions. 

Figure n shows that the excited state lifetimes depend strongly on the number of carriers 
(note the different vertical scale of each panel), the shorter lifetimes being shown by the 
N — 1 case. As a matter of fact, the N = 1 transition, which corresponds to the p —>■ s 
relaxation, represents the independent-particle limit of the ME cases shown in the same 
figure. For example, for N = 2 an independent particle filling gives a (M = 0,5 = 0) 
ground state, with the two electrons occupying the s spin-orbit als, and a (M = 1, S — 0) 
excited state, with one electron in the s orbital and another in the p orbital (see electronic 
configuration diagrams in Fig. I2^a)). Thus, the transition involves a one-electron scattering 
p — > s orbital, while the other electron remains as a spectator. Similar reasonings apply 
to iV = 3 — 5. Therefore, Fig. ^ shows that the Coulomb-interaction-free relaxation rate 
(N = 1 panel) generally gives an upper-bound to the actual rate when electron-electron 
interaction is taken into account. This result holds for all number of particles studied and 
most confinement strengths, although the trend is non-monotonic with N. 

The reduction of the relaxation rate noted above can be explained in terms of SE 
configurations mixing. In order to illustrate this, we analyze in detail the N = 2 and 
N = 3 cases in Fig. El where panel (a) represents the electronic configurations of the first 
excited and ground states which follow from an independent-particle filling, while panel 
(b) represents the two most important configurations contributing to the same states when 
Coulomb interaction is included. By comparison, one can see that in the non-interacting 
picture only the p — > s transition takes place, whereas in the interacting picture it is 
partially replaced by the d — > s and d — > p transitions. We then compare the relaxation 
rates of the individual SE scattering processes [panel (c)], fixing the transition energy 
Eg = hujQ in all cases in order to ensure that the comparison considers orbital effects only. 
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One can see that p — > s is the fastest transition. Therefore, when it is partially replaced 
by d — > s and d —>■ p the overall relaxation rate is reduced. Obviously, the stronger the 
mixing of configurations the larger the reduction. This explains why the N = 3 scattering 
rate is well below that of iV = 2: the ground state electrons already occupy the p-shell 
and the kinetic energy difference with respect to the excited states is then much smaller, 
which allows stronger Coulomb-induced mixing. Analogous logic can be used to explain the 
scattering reduction in the N > 3 systems. 

We next investigate the effect of a magnetic field B, applied along the vertical direction 
of the QD, on the ME relaxation rate. The magnetic field is expected to introduce new 
physics because it strongly modifies the SE energy levels, which now draw the well-known 
Fock-Darwin spectrum.— In particular, the states involved in the fundamental transition 
converge to the same (lowest) Landau level. This has important implications on the energy 
of the emitted phonon -which is reduced-, the regime of electronic correlations -which become 
stronger- and the SE configurations of the low-lying ME states -which differ from those at 
zero magnetic field-. In FigureEJwe plot the relaxation rate corresponding to iV = 2 in a QD 
with lateral confinement Hujq = 2 meV and width L z = 10 nm. We show the fundamental 
transition both in the singlet (S = 0) sector (solid lines) and in the triplet (S = 1) sector. 
We also compare the interacting and non- interacting case (thick and thin lines, respectively). 

It can be seen in the figure that the shape of the ME curves is again qualitatively similar 
to that of the SE ones. From previous investigations, we know that the scattering rate at 
fields exceeding a few Tesla is largely determined by the PZ scattering channel.— In addition, 
we note that the correlation-induced reduction of the relaxation rates changes with the 
field. Indeed, the factor of reduction tends to decrease as B increases, and at some point 
(B m 10.5 T for the singlet sector, B « 13 T for the triplet sector) the effect of correlations 
is reversed, so that a small enhancement is found instead of the more common reduction. 
This reversal in the behavior can be explained in terms of the mixing of configurations in the 
presence of a magnetic field. When no Coulomb interaction is considered only the p + — > s 
SE transition contributes to the total scattering, as shown in Fig. Ufa). However, if we take 
into account the Coulomb- induced mixing between the two most important configurations, 
Fig. E[b), new relaxation channels are opened, notably the d + —* s and d + — > p + SE 
transitions. The increasing relative weight of these transitions with B, along with the fact 
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that the d + — > p + transition becomes faster than the p + — > s one at sufficiently strong 
fields [Fig. life)], justify the observation that electronic correlations at high magnetic field 
yield slightly increased relaxation rates.— Yet, it may be worth pointing out that this takes 
place when the states involved in the transition have almost converged in energy (see inset 
in Figure EJ). Indeed, at such magnetic field values the investigated transitions are not the 
fundamental transition anymore, since higher-angular-momenta states have already come 
down in energy (this occurs at about B = 4 T in the QD we study, when the (M = 2, S = 1) 
level becomes the ground state). 

Recently, it has been suggested that electron-phonon coupling can be tailored in 
weakly-confined QDs and CQDs, so that SE lifetimes may be increased by orders of 
magnitude . 7 ! 11 ! 12 ! 13 The physical idea behind this possibility is to achieve an anti-phase rela- 
tion between the phonon wave and the electron wave function along the growth direction of 
the QD structure, i.e., to make the phonon wavelength along z a divisor of the quantum well 
width, thereby reducing strongly the value of the form factor in Equation El The phonon 
wavelength can be controlled through the lateral confinement, which is in turn determined 
either by electrostatic fields or by an axial magnetic field. The latter reduces the energy 
splitting between the low-lying SE states and thus increases the phonon wavelength in a 
controllable manner .''' 12 ' 13 The possibility of using external fields to suppress charge relax- 
ation also in ME systems may pose a significant support for the eventual fabrication of ME 
QD devices with small decoherece rates. However, the applicability is not straightforward: 
indeed, as mentioned before, the energy of the emitted phonon corresponding to the funda- 
mental transition decreases with increasing N, owing to the larger density of states.— As 
a result, the phonon wavelength may be too long to ever match anti-phase relation with 
the electron wave function in the growth direction. In order to explore this issue, in Fig. E] 
we represent the charge relaxation rate in a QD with Tiujq = 5 meV and L z = 15 nm with 
N = 2 and N = 3 electrons. For N = 2, two scattering minima show up at B = 0.6 T and 
B = 2.4 T. These are the same values as found in the N — 1 case, 7 which suggests that the 
N = 2 scattering is well described as an independent-particle event. However, the N = 3 
picture only shows one minimum at B — 0.55 T (the dip at B ~ 4 T is simply due to the 
crossing between the (M = 1, S = 1/2) and (M = 2, S = 1/2) energy levels). This is a clear 
manifestation of phonon energy reduction due to Coulomb interactions. For a larger number 
of electrons (not shown), when the density of states is still higher, the magnetic field brings 
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about frequent changes of energy levels with different symmetry, so that tuning the emitted 
phonon energy is no longer feasible. We then conclude that 5-induced suppression of charge 
relaxation rates holds for QDs with a small enough number of electrons, but it rapidly loses 
efficiency as the number of particles increases, owing to the increasing density of states. 

IV. COUPLED QUANTUM DOTS 

In this section we study charge relaxation rates corresponding to A-electron isospin tran- 
sitions in vertically CQDs as a function of the interdot barrier thickness The two QDs 
are supposed to be identical.— We focus on the transition between the lowest symmetric 
(G = 0) and antisymmetric (G = 1) solutions of the double quantum well with the spin 
quantum numbers of the ground state, i.e., the fundamental isospin (interdot) transition. 
This is indeed the fundamental transition of the system when the tunneling energy is smaller 
than the lateral confinement energy of the constituent QDs. 

Figure IH1 illustrates the results for a CQD structure with L z = 5 nm, huo = 5 meV and 
N = 1, 2, 3 electrons. Solid lines indicate the total scattering rate, while dashed and dotted 
lines represent the contribution coming separately from DP and PZ interactions, respectively. 
As for the intradot transition case, it is worth noting that the N = 1 calculation represents 
the independent-particle limit of the N = 2 and N = 3 systems. Since the maximum 
scattering rate of the N = 2 (N = 3) system is smaller than that of N = 1 by a factor of 
about 4 (2), we conclude that interdot transitions also benefit from a correlation-induced 
reduction of the scattering rates. Such reduction can be interpreted in terms of configuration 
mixing in analogous manner to the intradot case described in the previous section. 

We note in Fig. |H1 that the shape of the scattering rate curve is qualitatively similar for 
N = 1 and N = 3, with a dominating and oscillating DP scattering intensity at small barrier 
thickness, and a PZ contribution which becomes dominant as the barrier thickness exceeds 
~ 10 nm.i In contrast, the shape significantly differs for N = 2, where the PZ contribution is 
missing. Actually, a close inspection reveals that it has been suppressed by almost four orders 
of magnitude. This striking result cannot be interpreted as a consequence of correlations on 
the electron wave function, as in previous sections, because the orbital part of the electron- 
phonon interaction matrix element is identical for DP and PZ scattering mechanisms (see 
Eq. OJ), and therefore the effect should be apparent also for DP scattering. Therefore, the 



10 



origin must be connected with the effect of Coulomb interaction on the emitted phonon 
energy. 

To better understand this result, in Fig.[7|we compare the L b dependence of the tunneling 
energy, and hence that of the emitted phonon, for N = 2 with and without Coulomb 
interaction. Clearly, Coulomb interaction is responsible for a significantly faster reduction 
of the tunneling energy. Consequently, the values of Lb which would lead to maximum 
PZ scattering for non- interacting electrons (e.g., L& ~ 11 nm), are now associated with 
very small phonon energies. As a result, the phonon density of states is very small and the 
electron-phonon coupling is strongly reduced. This behavior is no longer found in the N = 3 
system. An interpretation of this observation is provided in the diagrams of Fig. where 
we show the dominant SE configurations for the interacting and non-interacting electrons. 
In the absence of Coulomb interaction, the only configuration for the G = state has two 
electrons in the lowest SE symmetric orbital (cr 2 ). Similarly, for the G = 1 state the only 
significant configuration is that with one electron in the lowest SE symmetric orbital and 
another in the lowest antisymmetric orbital {a a*). However, when Coulomb interaction is 
included for G = we mostly obtain a linear combination of a 2 and (cr*) 2 configurations. The 
thicker the barrier, the smaller the tunneling energy and the larger the mixing between these 
two configurations. In the limit where both configurations have equal weight, the G = 
state will be degenerate with the G = 1 one. In other words, Coulomb interaction tends to 
make G = and G = 1 solutions converge in energy with increasing Lb- This is not possible 
for N = 3 because the unpaired electron prevents a similar mixing of configurations which 
use a and cr* orbitals while preserving the total parity quantum number. Therefore, the only 
significant electronic configuration for G = and G = 1 are cr 2 er* and a(a*) 2 , respectively. 
The transition between these two configurations can be envisaged as the relaxation of a 
single hole from cr* to cr orbitals, which explains the similar behavior of N = 3 as compared 
to iV = 1. 

V. SUMMARY 

We have investigated the effect of Coulomb interaction on the charge relaxation rates 
of disk-shaped QD systems with few electrons, due to incoherent coupling with acoustic 
phonons. We have studied both intradot transitions in single QDs and interdot (isospin) 
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transitions in vertically CQDs. Coulomb interaction affects the scattering rates in two ways: 
first, it changes the emitted phonon energy and momentum, and second, it changes the elec- 
tron wave function through electronic correlations. Both effects have a significant influence 
on the electron relaxation rate in a way that generally leads to a reduction of the intradot 
and interdot transition rates. This trend is gradually quenched by axial magnetic fields. On 
the other hand, the increasing density of states with higher number of particles reduces the 
energy of the emitted phonons, which renders PZ interaction increasingly important, so that 
it may eventually become the dominant scattering mechanism even at zero magnetic field. 
We have finally shown that the suppression of electron scattering in weakly- confined QDs, 
which was recently suggested for SE systems,— ilZill also applies to multi-electron QDs with 
a sufficiently small number of particles. 
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Figure U because it occurs for characteristic frequencies hu>o < 0.1 meV, but it is apparent e.g. 
in Figure Hfc). 

24 The large shift from N = 2 to N = 3 is due to the occupation of the p-shell by the ground state 
electrons, which greatly reduces the kinetic energy difference with respect to the excited states. 

25 It may be worth pointing out that the phonon energy of the plotted transitions in Figure 0] is 
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identical, so the different rates are due to orbital effects only. 

Note that the role of electronic correlations is secondary for this discussion since we are con- 
cerned about the electron wave function in the vertical direction only, which is the first solution 
of the quantum well for all low- lying QD states due to the strong vertical confinement. 
The scattering rates of SE CQDs with slightly anisotropic QDs have been studied in Ref. 
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FIG. 1: Charge relaxation rate vs lateral confinement for a QD with vertical width L z = 10 
nm filled with N = 1 to N = 5 electrons. Solid lines: total scattering rate. Dashed line: DP 
contribution. Dotted lines: PZ contribution. Note the different vertical scale in each panel. Next 
to the right axis, we indicate the quantum numbers (M, S) of the states involved in the transition. 
The downward arrows in each panel point at the confinement energy leading to emitted phonon 
energies of 1.3 (solid arrowheads) and 2 meV (empty arrowheads). Inset in the top panel: SE 
configurations involved in the N = 1 transition. 
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FIG. 2: SE electronic configurations (Slater determinants) involved in the fundamental transition 
for the QDs of Fig. ^ for N = 2 and N = 3 for (a) non-interacting and (b) interacting electrons. 
In the interacting case only the two most weighted configurations are shown, (c) SE scattering 
rates vs lateral confinement energy between selected Fock-Darwin orbitals at fixed phonon energy 
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Fig.3 
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FIG. 3: Charge relaxation rate vs magnetic field for the lowest singlet (solid lines) and triplet 
(dashed lines) transitions for N = 2 in a QD with L z = 10 nm and hujo = 2 meV. Thick lines: 
interacting case. Thin lines: non-interacting case. Inset: lowest-lying energy levels and their 
quantum numbers (M,S); arrows indicate the transitions under study. 
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FIG. 4: SE electronic configurations (Slater determinants) involved in the fundamental transition 
for N = 2 and N = 3 for (a) non-interacting and (b) interacting electrons in a magnetic field. In 
the interacting case only the two most weighted configurations are shown, (c) SE scattering rates 
between selected Fock-Darwin orbitals vs magnetic field. QD parameters are the same as in Fig. 01 
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Fig.5 
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FIG. 5: Charge relaxation rate vs magnetic field for N = 2 and N = 3 in a QD with L z = 15 nm 
and hujQ = 5 meV. Insets: lowest-lying energy levels along with their quantum numbers (M,S); 
arrows indicate the transition under study. 
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Fig. 6 
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FIG. 6: Charge relaxation rate vs interdot barrier thickness in a CQD with L z = 5 nm, ftwo = 5 
meV and iV = 1,2,3 electrons. Solid lines: total scattering rate. Dashed line: DP contribution. 
Dotted lines: PZ contribution. Note the different vertical scale in each panel. Next to the right 
axis, we indicate the quantum numbers (M, G, S) of the states involved in the transition. For 
N = 2, the PZ contribution cannot be distinguished on this scale. Insets: lowest-lying symmetric 
(solid line) and antisymmetric (dashed line) energy levels. 
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FIG. 7: Energy of the emitted phonon vs interdot barrier thickness for N = 2 interacting (solid 
line) and non-interacting (dashed line) electrons in the CQD structure of Fig. H3 The diagrams 
illustrate the dominant SE configurations with and without Coulomb interaction. 



20 



